Quantum criticality of a one-dimensional Bose-Fermi mixture 
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The one-dimensional interacting Bose-Fermi mixtures, exhibiting quantum phase transitions at 
zero temperature, are particularly valuable for the study of quantum critical phenomena. In the 
present paper, we analytically study quantum phase diagram, equation of state and quantum criti- 
cality of the Bose-Fermi mixture using the thermodynamic Bethe ansatz equations. We show that 
thermodynamical properties display universal scaling behaviour at quantum criticality. Further- 
more, quantum criticality of the Bose-Fermi mixture in an harmonic trap is also studied within 
the local density approximation. We thus demonstrate that the phase diagram and critical proper- 
ties of the bulk system provide insights into understanding universal features of many-body critical 
phenomena. 

PACS numbers: 03.75.Mn, 03.75.Hh, 02.30.Ik, 05.30.Rt 



I. INTRODUNCTION 

Mixtures of ultracold bosonic and fermionic atoms 
have attracted intensive studies both experimentally and 
theoretically |3-[3:]. By loading cold atoms in one- 
dimensional (ID) waveguides and tuning the effective in- 
teractions by Feshbach resonance, it is possible to sim- 
ulate striking quantum many-body phenomena in ID 
strongly correlated systems in the whole regime of in- 
teraction strength [3t9- The exquisite tunability with 
ultracold atoms confined to low dimensions has provided 
unprecedented opportunities for investigating and testing 
the theory of exactly solvable many-body systems (7l-[TT| . 
These include remarkable experimental progresses in the 
realization of Tonks- Girardeau gas [Ti, iSj], super- Tonks- 
Girardeau gas 0] , Yang- Yang thermodynamics for ultra- 
cold Bose gas of ^''Rb [10], and the exotic density profiles 
of the attractive Fermi gas in an harmonic trap JJj] . The 
current experimental progresses are capable of simulating 
ID Bose-Fermi mixtures. 

Recently, various theoretical methods have been used 
to study quantum phases and correlations of the ID Bose- 
Fermi mixtures, such as the mean-field approach 11211 
and Tomonaga-Luttinger liquid (TLL) theory ^isl. Il4| . 
etc. The ID Bose-Fermi mixture with equal masses of 
bosons and fermions and with the same strength of delta- 
function interaction between boson-boson and boson- 
fermion is exactly solvable |15| . This model has been 
attracted a renewed interest |16l - [l9| due to the experi- 
mental progress with cold atomic systems. 

Particular theoretical interest has been paid to the 
ground-state properties at zero temperature 
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However, there are very few studies on the thermodynam- 
ics and quantum critical phenomena of the model. In re- 
view of the realistic cold atomic systems trapped in exter- 
nal potentials at finite temperatures, it is significantly im- 
portant to understand how to unambiguously determine 
the zero temperature phase diagram from the knowledge 
of finite temperature quantities of trapped gases. In the 
ID mixture of quantum gases, true quantum phase tran- 
sitions occur as the driving parameters vary across the 
phase boundaries at zero temperature, such as chemical 
potential, magnetic field and densities, etc. In particular, 
ID quantum critical phenomena associated with quan- 
tum phase transitions at zero temperature give physical 
origin of quantum criticality (20l . [2l| and provide an in- 
sight into understanding of universal scaling behaviour of 
thermodynamical properties in quantum critical regimes 
(22h25i l. By using the universal scaling functions, it has 
been demonstrated that the zero-temperature phase di- 
agrams of various systems can be mapped out from the 
finite-temperature density profiles [22, 24, ^fij- Most re- 
cently, the high-resolution imaging techniques have al- 
lowed to measure the density profiles and density fluctu- 
ations of the trapped atomic gases very precisely [l^-dll 
and thus provide essential tools to study quantum phase 
transitions and quantum criticality. 

In general, quantum fluctuations are strongly cou- 
pled with thermal fluctuations in the quantum critical 
regime. Therefore, quantum criticality is among the 
most challenging problems in condensed matter physics. 
In order to extract correct universal scaling functions 
which control proper thermal and quantum fluctuations 
at quantum criticality, high precision of flnite temper- 
ature thermodynamics is desirable. Usually, accessing 
to the thermodynamic properties of intcgrable models at 
finite temperatures is notoriously difficult and presents 
a formidable challenge in theoretical and mathematical 
physics. In the present paper, we analytically determine 



2 



the zero temperature phase diagram of the integrable 
Bose-Fermi mixture. We further derive equation of state 
and explore universal scaling behaviour of thermodynam- 
ical properties at quantum criticality using the thermo- 
dynamical Bethe ansatz (TBA) equations. Using exact 
analytical result obtained, we also demonstrate that the 
zero temperature phase diagram and quantum critical- 
ity can be mapped out from finite temperature density 
profiles of the trapped gas within the local density ap- 
proximation. 

The paper is organized as follows. In section we 
present the TBA equations for the model and analyt- 
ically determine the phase diagram of the Bose-Fermi 
mixture at zero temperature. In section lllli we derive 
equation of state and explore the universal scaling be- 
havior of the density and compressibility near the critical 
points. In section ITVl the quantum criticality of the gas 
in an harmonic trap is studied within the local density 
approximation. A conclusion is presented in the section 
lYl 



II. MODEL AND PHASE DIAGRAM AT ZERO 
TEMPERATURE 

We consider a ID interacting Bosc-Fcrmi mixture de- 
scribed by the Hamiltonian 



H ^ I dx . 

Iq \2mi, 



" 2m f J 



(1) 

where ^b, ^/ are boson and fermion field operators, 
nib, TTif are the masses, /if,, /i/ are chemical poten- 
tials of bosons and fermions, and gbb, gbf are boson- 
boson and boson-fermion interaction strengths, respec- 
tively. Here we consider the fully-polarized fermions, 
therefore, the Pauli principle excludes their s-wave in- 
teraction {gff = 0). This model is exactly solvable 
[l6| for equal masses and equal boson-boson and 
boson-fermion interaction strengths, i.e., rub = mj = m 
and gbb = gbf = g- For our convenience, we assume 
h = 2m = 1. The first quantization form of the exactly 
solvable Hamiltonian ([IJ can be written as 



N 



H 



^ dx. 



x,)-t^N~^iNf-Nb) 
(2) 

with c = rng/ti^. Here the particle number N = Nb + Nf 
with Nb bosons and Nf fermions. The chemical poten- 
tial /i and the effective magnetic field H are defined as 
fi — (/i/ -I- fib) /2 and H = fj,f — jib- The many-body 
wave function is supposed to be symmetric with respect 
to indices i = {l,2,...,A^f,} (bosons) and antisymmetric 
with respect to {Nb + l,Nb + 2, TV/} (fermions). Thus 



the TV-body wave function can be written as N\ x N\ su- 
perpositions of individual plane waves associated with N 
quasi-momenta kj with i = 1,...,A^ by means of the 
Bethe ansatz 



15|,ll6|. 



The spectrum of the system is given hy E = '}2^=i 
where the quasi-momenta ki is subject to the so-called 
Bethe ansatz equations (BAEs) [H, [ll] . In thermody- 
namic limits, i.e. N ^ oo, L ^ oo and N/ L is finite, and 
at finite temperatures, the equilibrium states become de- 
generate. Yang and Yang [S^l showed that true physical 
states in integrable systems can be determined from the 
minimization conditions of Gibbs free energy subject to 
the Bethe ansatz equations. Accordingly, in the thermo- 
dynamics limit with Nb/ L and Nf/L fixed, minimization 
of the Gibbs free energy gives rise to the following non- 
linear integral equations [l9|, i.e., the TBA equations for 
the integrable Bose-Fermi mixture: (JtB — 1) 

e (k) = k^ — Hf 

/oo 
ai (A - k) In (1 + exp {-ip (A) /T)) dA, 
-OG 

V (A) = fif- Hb 

/>OG 

-T I ai (fc - A) In (1 + exp (-e (k) /T)) dk, (3) 

For fixed temperature T 



1 ic 

2tt (icy-'/i+x-' • 

and chemical potential /i/, /if,, the pressure is given by 



where ai {x) 



T 



In 1 + exp 



e(fc) 
T 



dk. 



(4) 



The particle density and compressibility for fermions and 
bosons, and the entropy per length can be obtained from 



S/L 



N,/L = dp/ dm 

drii/dfi, 
dp/dT. 



f,b 



(5) 
(6) 
(7) 



The phase diagram of the Bose-Fermi mixture can be 
analytically determined from the TBA equations (jS]) in 
the zero temperature limit. As shown in Fig. [TJ the 
phase diagram consists of three quantum phases: pure 
bosons, pure fermions and the mixture of bosons and 
fermions except the vacuum, which are separated by four 
boundary lines with condition nf = or Ub — 0. In 
order to obtain the full phase diagram of the Bose-Fermi 
mixture, we introduce two sets of the TBA equations 
with different reference states [s^. The TBA equations 
based on the fermionic reference state determine the two 
boundary lines for H > 0, whereas the TBA equations 
based on the bosonic reference state determine the two 
boundary lines for H < 0. 

In the limit T — )■ 0, the TBA equations ^ based on 
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FIG. 1: (Color online) Phase diagram in fi — H plane. Three 
distinguished phases are resulted from varying the chemical 
potential and magnetic field, i.e. pure boson phase for H < 
and jj, > H/2; pure fermion phase below the phase boundary 
(|12|l in the region fj, > —H/2; and the mixture of bosons and 
fermions above the phase boundary p2|l in the region H > 0. 



the fermionic reference state reduce to 

fc2 



e{k) 



ai {A~ k)ip (A) dA, 



ip{A) = Hf - ^ib+ ai(fc-A)e {k)dk, (8) 

where the dressed energies e~{k) and (f~{K) correspond 
to the occupied states for k € [— fcF,fc_F] and A e 
[—Kp, Ap] , respectively. The positive parts of the dressed 
energies e{k) and <p(A) corresponds to the unoccupied 
states. The integration boundaries kp and Ap charac- 
terize the Fermi points with the conditions (ikp) — 
<y9~ (±Ai?) = 0. The pressure is given by 



P - 



1 

'2^ 



{k)dk 



(9) 



at zero temperature. We calculate the particle densities 
through the relations ([5]). Nevertheless, the boundary 
lines, which correspond to Uf = or rib = 0, can be 
determined by analyzing the dressed energy at the point 
kp = OT Ap = 0. 

For a pure fermion state, ip (A) is gapfull, i.e. ip (A) > 
0. Thus the TBA equations (HI reduce to the free fermion 
potential e{k) = k^ — jif. The phase boundaries for 
the phase transitions from vacuum into the pure fermion 
state and from pure fermion state into the mixture of 
bosons and fermions are determined by 



e"(0) < 0, 
ip~ (0) < 



(10) 

(11) 



respectively. From the condition (|10p . we obtain the 
phase boundary for pure fermions /i/ > 0, or equivalently 
/i > -H/2. 



From the conditions (jlip and e {k) 
the phase boundary 



- 2tt 



{Afif + 1) arctan - 



fc^ — /z/, we obtain 



(12) 



for the coexistence of bosons and fermions. Here we used 
the dimensionless units H = H/eq and flf = fif/so with 
eo = c^- For strong coupling regime, i.e. fj,f <^ Eq, or 
H So, find phase boundary condition (|T2|) reduces 
to 



> 



2/3 



2 / ~\2/3 

15 (^-^ 



2H 



(13) 



with jl — fJ./eo. For weak coupling regime, i.e. fif Eq, 
or H ^ Eq, the phase boundary condition becomes 



H 



2 

TT 



(14) 



On the other hand, at zero temperature, the TBA 
equations with the Bose state as the reference state is 
given by [s^ 



e{k) = k^-Hb+ / a2{k- k')e- {k')dk' 

J — kp 




+ Ui (fc - A) (A) dA 



if (A) = fih-^if- ai(fc-A)e {k) dk (15) 

J-kp 

with the Fermi points e~ (ikp) = (p^ {±Ap) = 0. Simi- 
larly, the dressed energies e~ (fc) and ip^ (A) correspond to 
the occupied states for k G [—kp, kp] and A e [—Ap, Ap], 
respectively. We see that for the Bose reference state the 
fully-polarized fermions provide a ferromagnetic order- 
ing at the ground state. If < 0, i.e. fif < fit, the 
dressed energy ip{A) is greater than zero. The dressed 
energy (p is gapful. Thus the ground state is a pure bo- 
son state. Therefore, for H < 0, the TBA equations (fT5)) 
reduce to Yang- Yang thermodynamics equations for the 
Lieb-Liniger Bose gas [32| 

fkp 

e{k) ^ k^ - fih+ / 02 {k - k') €- {k') dk' 



from which we easily determine the phase boundary for 
the phase transition from vacuum into the pure boson 
state, i.e., ^ > H/2, see Figure[TJ 



III. EQUATION OF STATE AND UNIVERSAL 
SCALINGS 



Recent experiments on quantum criticality of ultra- 
cold atoms jH, [3l| and theoretical scheme of mapping 
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out quantum criticality of cold atoms [22| - |2a | open the 
possibility to explore such universal behavior in low di- 
mensional multi-component interacting Fermi and Bose 
gases. As the temperature is tuned over the crossover 
temperatures the scaling functions of thermodynamical 
properties give rise to universal behavior which entirely 
depends on the symmetry of excitation spectrum and di- 
mensionality of the system. This gives a promising way to 
explore the hidden symmetry of these models, for exam- 
ple, the quantum Ising model with transverse file displays 
emergent E8 symmetry [35|. In the critical regime, the 



thermodynamic functions of the homogeneo us g as can be 
cast into some universal scaling forms [1^ l2l| For ex- 
ample, the density and compressibility near the critical 
point II = fic can be written as 



n {T,iji) ^ no + T 



-+1-7 



(16) 



(17) 



Here uq (ko )is the regular part of the density (com- 
pressibility) induced from the background, Q (J^) is a 
universal scaling function describing the singular part of 
the density (compressibility) near the critical point /ic, d 
is the dimensionality of the system, z is the dynamical 
critical exponent and v is the correlation length expo- 
nent. From the above relation, the dimensionless univer- 
sal scaling functions G ( '^^^Jt" ) and J^{^^f^) display uni- 
versal scaling behaviour near the critical point of /i = /i^, 
i.e. the density(compressibility) curves with a subtrac- 
tion of the background density (compressibility) intersect 
at the critical point for different temperatures. This fea- 
ture can be used to detect the phase boundaries at zero 
temperature from finite temperature density profiles of 
the trapped gas. 

Before discussing the universal scaling behavior of 
the Bose-Fermi mixture, we will discuss quantum crit- 
icality of several simple examples. The simplest ex- 
ample is the ID ideal Fermi gas. The ideal Fermi 
gas obeys the Fermi-Dirac distribution, it is easy 
to derive the density distribution for free fermions 
nfiT,fj.) = -Vr/(2v^)LH/2 (-exp(^)), where 
Lin {x) = ^'iZi x'- /l^ is the standard polylogarithm func- 
tion. In comparison with the universal scaling Ea.(|16p. 
the critical exponent z — 2 and the correlation length 
exponent v — 1/2 with the dimensionality d = \ can be 
read off the universal scaling form. Here the scaling func- 
tion is Q (x) = —1/ (2-^7?) Lii/2 {— exp (x)). There is no 
background density for the vacuum-Fermi-gas transition, 
i.e. no = 0, where the critical point /ic = 0. It was also 
shown that d = 1 , z = 2 and = 1 /2 for the vacuum into 
the TLL phase transition in the ID hard-core bosons [2^ . 
the ID attractive Fermi gas [l^] and the ID interacting 
Bose gas with strongly repulsive interactions [1^. 

In order to investigate the quantum critical behavior 
in the vicinity of phase boundary between the phase of 
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FIG. 2: (Color online) Contour plot of entropy S vs chem- 
ical potential from the TBA ((3|: quantum criticality driven 
by chemical potential for H — O.leo- The crossover tempera- 
tures (white squares and triangles) separating Vacuum, TLLf 
and TLLm from the quantum critical regimes are determined 
from the breakdown of linear temperature-dependent entropy 
from (|22|) . Here the TLLf stands for the TLL of fermions 
where exponentially small number of bosons are populated at 
finite temperatures. Whereas TTLm denotes the TLL of the 
mixture. The vacuum evolves into a quasi-classical regime at 
finite temperature. 



a mixture of bosons and fermions and the phase of the 
fully-polarized fermions, we will derive the equation of 
state from the TBA equations ([3]). For strong interacting 
regime, i.e. c ^ 1, or T/eo ^ 1, we can rewrite the TBA 
equations ^ 



e{k) 
<^(A) 

where 



k^/fi - A 
M/ - Mb - 



(18) 



Ape 



4A2 



T 



2Tc 



P = 1 



P2 



I3T' 



Lis 



ai (A) In ( 1 

^ 4c2-48A2 

-oo (c2 + 4A2)3 

A 



4c (4c2 - 48A2) p2 
(c2 -I- 4A2)^ ' 

VP (A) 



exp 



■In 1 



T 



dA, 



dA, 



(19) 



With the help of these relations, the pressure (|4]) can be 
calculated in a straightforward way 



V 47r 2 V / 



(20) 



that serves as the equation of state of the model with 
strong repulsion. The thermodynamical properties can 
be obtained from the usual thermodynamical relations. 
This analytical equation of state ([20]) essentially cover the 
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universal TLL thermodynamics and encode the critical 
exponents in the critical regimes. 

At very low temperatures, i.e. T ^ Sq, the thermo- 
dynamics of the model is governed by the TLL physics 
associated with a linear dispersion. In the mixed phase 
of bosons and fermions, the low energy physics belongs 
to a universality class of a two-component TLL [itI [18| . 
In this low temperature limit, we further calculate the 
pressure 



2AV2 



8 



where 



A 



H A fH p 
2HAo cttT^ 



tan 



2An 



TTC 



12HA, 







(21) 



with Aq = cy He ~ 1- After a lengthy algebra, we find a 
universal leading order of temperature corrections to the 
free energy 



ttCT^ ( 1 1 



6 Vj 



(22) 



where the two velocities in strong repulsive regime are 
given by 



■sin(7ra), 
( 4 

Vf — 27rn ( 1 (ttq; -f sin(7ra)) 



(23) 



The parameter a is determined by the relation a sa nhjn 
for a small H <^ 1. The ground state energy is given by 



7 V2 



sin(7rQ:) 



a) cos(7rQ;] 



(24) 

The TLL is maintained under a crossover tempera- 
ture at which the linear temperature-dependent entropy 
breaks down, see Fig. [2j The exact analytic expression 
of thermodynamic functions ([2^ from the usual TLL de- 
scription is only accurate for a limited range of temper- 
atures and density. However, the TLL description is in- 
capable of describing quantum criticality for it does not 
contain the right fluctuations in the critical regime. The 
equation of state (j20l) contains proper universal scaling 
functions which control full thermodynamical properties 
in the quantum critical regimes. Near critical point, the 
thermal dynamical properties can be cased into universal 
scaling forms, for example, (|16p and p7)) . There exists no 
longer free fermions in the TLLp phase in Fig l2] due to 
the presence of exponentially small number of bosons at 
finite temperatures. They behave like a TTL for temper- 
atures below the crossover temperatures (white triangles 
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FIG. 3: (Color online) The density ni, and compressibility 
vs chemical potential ^ for H/eo = — 1 at different tempera- 
tures. The curves intersect at the critical point /i = H/2, i.e. 
the phase boundary between the pure boson phase and the 
vacuum. Here t = T/eq. 



and squares). The TLL^i stands for a two-component 
TLL of the mixture of bosons and fermions described by 

It is straightforward to work out the critical exponents 
for the phase transition from vacuum into the the free 
fermions and into the Bose gas phase i.e. d = 1, z = 2 
and = 1/2, see [ll,!!!!!^. Using the TEA equations 
([3]), we verify the scaling behavior of the mixture for the 
phase transitions from vacuum into the pure Bose state 
in figure [3] and from vacuum into the pure Fermi state 
in figure |4l We see that both the density and compress- 
ibility at different temperatures intersect at the critical 
point. For both cases, the compressibility always evolves 
a round peak as temperature decreases. This indicates 
the density of state changes when the phase transition 
occurs. As the temperature approaches the limit T — > 0, 
the compressibility tends to divergence. 

Comparing the density, which is numerically obtained 
from the pressure ([4]), with the universal scaling form of 
Ea. (|16p . we can extract the critical exponent z = 2 and 
the correlation length exponent v — 1/2 with the dimen- 
sionality d = 1. Here we numerically demonstrate the 
universal scaling behavior of the mixture of bosons and 
fermions in Figs [31 S] and O For practical convenience, 
here we have chosen the density of bosons rih to demon- 
strate the intersections for the phase transitions from vac- 
uum into the phase of pure bosons and from the phase 
of fermions into the mixture, where ni, does not have a 
background near the transition points, see Figs. [3l and 
[5] For a fixed effective magnetic field H, we see that by a 
proper temperature scaling the density curves nb(T)/y/T 
at different temperatures intersect at the points /ic. In 
the Fig. [5l we display the scaled density distributions for 
different temperatures by numerically solving the TBA 
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FIG. 4: (Color online) The density nf and compressibility Kf 
vs chemical potential jj, for H/eo = 1 at different tempera- 
tures. The curves intersect at the critical point p — —H/2, 
i.e. the phase boundary between the pure fermion phase and 
the vacuum. Here t = T/eo- 
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FIG. 5: (Color online) The density nt and compressibility 
vs chemical potential fi for H/eo = 1 at different tempera- 
tures. In this setting, the intersection nature can map out 
the zero temperature phase boundary for the phase transi- 
tion from the free fermions into the mixture of bosons and 
fermions. Here t — T/eo- 



equations (I3|)-(l5]). It is also seen that the compressibility 
curves at different temperatures intersect at the critical 
point ^c- The compressibility tends to divergent as the 
temperature tends to zero. It evolves to a round peak 
at low temperature due to the change of the density of 
states around the critical points. The quantum criticality 
near the critical points associating the phase transitions 
from vacuum into the phase of pure bosons and the phase 
of pure fermions reveals a subtle difference in thermody- 
namical properties, see Figs |3] and [S) 
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FIG. 6: (Color online) The density distribution at different 
temperatures with N^a\j^/a? = 1 and a — 0.5 . 














1=0.002 








t=0.004 






* ^ , N 


1=0.0075 






^"^""""""^---^ \ 

^""""""--•^ " ^ \ 




(b) 







1.20 

x/a,„N" 



FIG. 7: (Color online) The densities nt, and n/ vs normal- 
ized position for Na\j^/a^ = 1 and a = 0.5 at different tem- 
peratures. The density curves intersect at the critical point 
that maps out the phase boundaries of the trapped gas. Here 
t = T/eo. 



IV. QUANTUM CRITICALITY IN THE 
HARMONIC TRAP 



In experiment vifith cold atoms, the ID quantum gas is 
realized by tightly confining the atomic cloud in two (ra- 
dial) dimensions and weakly confining it along the axial 
direction in an external harmonic trap. For the mixture 
of bosons and fermions in a harmonic trap, we can cal- 
culate its density distribution profiles by evaluating the 
thermodynamical dynamics within the local density ap- 
proximation. According to the local density approxima- 
tion, the system reaches local equilibrium in each small 
intervals around each point x in the external trap. The 
density distribution of the trapped gas is then obtained 
via the local equation of state [l^ Is^l- Within the lo- 



7 



cal density approximation, the chemical potentials in the 
equation of state ^ as well as in the TBA equations are 
replaced by the local chemical potentials given by 

Hb{x) = fib{0)-Vb{x), (25) 
fif ix) = fif (0) - Vf ix) . (26) 

Here the external potential is defined as VJ, (x) = Vf (x) — 
muP'x^ jl with harmonic frequency a; and the character- 
istic length for the harmonic trap is a = a/ h/mui. For 
this case, equations ((25|) and ((26| can be alternatively 
represented as 

^^ (y) /eo = (0) /eo - 

for a fixed i7, where the dimensionless coordinate is given 
by y = x/{a'^c). From the Bethe ansatz equations, the 
dimensionless density nb/c and nj/c can be obtain for 
fixed dimensionless chemical potential fJ-b/so ^.nd Hf/sQ. 
The total particle number N is obtained from 

Nafj, _ ^ r Ub {y)^ns jy) ^ 



with the ID scattering length am — — 2/c. We define a 
polarization rate between the Bose particle number and 
the total particle number a = Nb/N. For fixed /i (0) /eq 
and H/eq, we can calculate Na\jj/a? and a. In the pres- 
ence of the confined potential, the length scale of the 
system at quantum criticality is still much smaller than 
the trap size. Therefore, the critical behavior of the ho- 
mogeneous gas can be mapped out by the density profiles 
of gas at finite temperatures [13] ■ 

We fix Na\jj/a? and a in the trapped gas of the mix- 
ture, the density profiles reveal a universal scaling behav- 
ior of the gas. Fig. [6]shows the density profiles of bosons 
and fermions in the harmonic trap for Na\jj/a? = 1 
and a = Nb/N = 0.5 at different temperatures. Here 
we find that bosons and fermions coexist in the trap 
center companied by the phase of pure fermions at the 
edges. We further demonstrate how to map out the zero- 
temperature phase boundaries from the density profiles 



of the trapped gas at finite temperatures. In Fig. [Tl^a) 
and (b) , we demonstrate the scaled density distributions 
of bosons and fermions. It is clearly seen that the scaled 
density curves for different temperatures intersect at the 
critical point of the trapped gas. Thus the critical point, 
separating the mixture of bosons and fermions from the 
phase of pure fermions, is mapped out. Similarly, the 
density curves of fermions intersect at the critical point 
that maps out the phase boundary for the phase transi- 
tion from vacuum into the phase of free fermions. 



V. CONCLUSION 

In summary, we have studied the phase diagram, uni- 
versal TLL and quantum criticality of the ID Bose-Fermi 
mixture by means of the TBA equations. We have de- 
rived the equation of state and universal TLL thermody- 
namics of the model for strong repulsion. We have proved 
that the low energy physics of the Bose-Fermi mixture 
are described by a two-component TTL. Universal scal- 
ing behavior of thermodynamical properties at quantum 
criticality provides a physical origin of quantum critical 
phenomena. Furthermore, the quantum criticality of the 
Bose-Fermi mixture in an harmonic trap has been stud- 
ied within the local density approximation. It turns out 
that the phase diagram and critical properties of the bulk 
system can be mapped out from the density profiles of 
the trapped mixture gas at finite temperatures. Our ex- 
act results can help with experimental study of quantum 
critical phenomena in a ID harmonic trap. 
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